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1. Introduction 

The study of possible modifications of hadron properties in the nuclear medium is one 
of the central problems of contemporary nuclear physics. In principle, these and related 
phenomena in nuclear physics are governed by the fundamental theory of the strong interac- 
tions, quantum chromodynamics (QCD). However, although QCD has been very successful 
in explaining a large class of hadronic processes at high energy and large momentum transfer, 
typical nuclear phenomena at lower energies cannot be derived from QCD with the theoreti- 
cal tools presently available. The difficulty of using QCD for phenomena at the nuclear scale 
is related to the nonperturbative nature of these. Due to the asymptotic freedom property 
of QCD, high energy processes are calculable by perturbative techniques in the quark-gluon 
coupling constant. On the other hand, since there are no reliable systematic approxima- 
tion schemes in field theory for performing nonperturbative calculations, the construction of 
models is an important aspect of low energy QCD. While there is considerable optimism that 
eventually one will be able to solve QCD numerically on the lattice using supercomputers, 
the development of analytical approximation methods are in urgent need to make contact 
with the wealth of data on nonperturbative phenomena presently available, or that will be 
available when the new experimental facilities under construction start operating. The S 
expansion is an example of a method recently developed aiming to study nonperturbative 
phenomena in field theory. 

The idea of the 6 expansion is to perturb the original theory by the introduction of an 
artificial expansion parameter 6, absent in the original theory. The parameter 6 is introduced 
in such a way that it interpolates between the theory one wants to solve and another theory 
that one knows how to solve. The 6 expansion can be formulated in two different forms, the 
logarithmic 6 expansion IQ and the linear 6 expansion 0- Q . In this paper we consider the 
linear form. Specifically, let C be the Lagrangian density of the theory one wants to solve, 
and Co the Lagrangian density of the soluble theory. Then, the interpolating Lagrangian 
density C{6) is defined as 

C{S) = (1 - S)Co + 5C = Co + 5{C-Co), (1) 

so that C{0) = Co, C{1) = C and Cq is a function of an arbitrary mass parameter /i. The 
next step involves the evaluation of desired physical quantities as a perturbation series in 
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powers of S, and then S is set equal to 1 at the end. A crucial aspect of the method is 
the recognition that Cq involves arbitrary unknown (dimensionful and/or dimensionless) 
parameters. 

Fixing the arbitrary parameter /i is the step which brings all nonperturbative information 
contained in the perturbative calculation. Several ways to fix the arbitrary mass parameter 
have been proposed in both versions of the 6 expansion as well as in the related methods. 
One physically appealing way to fix the unknown parameters, which is the one adopted 
here, is the principle of minimal sensitivity (PMS) introduced in Ref. 0. This variational 
principle amounts to the requirement that a physical quantity P(/i) should be at least locally 
independent of these parameters, which implies that 

^1.^0, (2) 

at (5 = 1. The solution to the PMS equation gives /i as a function of the original parameters 
of the theory including the coupling constant. The 6 expansion, together with the criterion 
of the PMS of physical observables, is known as the optimized 6 expansion. The convergence 
of the optimized 6 expansion has been proved in Ref. [0 for a quantum mechanical problem. 
The different forms of the 6 expansion have been successfully applied to many different 



problems in quantum mechanics 0, particle theory 0- ||T0[, statistical physics ||TT| and 
lattice field theory |^- |12|. Most applications show that one is always able to reproduce 
traditional non perturbative results already at lowest order in 6. 

Obviously different approximations give different prescriptions so as to select a subset 
of Feynman diagrams among the infinite set which describes a physical process and within 
the 6 expansion this selection is done in an essentially perturbative way. Also, as we shall 
explicitly see, the same order in 6 can contain diagrams which would belong to different 
orders if we were using other approximations. A drawback of traditional nonperturbative 
analitycal approximations is that one has to sum an infinite subset of graphs so as to consider 
all orders in the coupling. This procedure generates problems related to the inclusion of 
higher orders or nonperturbative renormalization or both. 

One advantage of the method presented here is that one deals with a reduced number 
of Feynman graphs so that renormalization can be carried out in a perturbative way before 
the PMS produces the final finite nonperturbative results. Also, because there is no self 



consistency involved, it can be considered more economical as far as numerical computa- 
tions are concerned. Motivated by these advantages, we have recently implemented the 



optimized 6 expansion for the Walecka model [Q. We have investigated vacuum effects by 
neglecting exchange diagrams and have shown that the relativistic Hartree approximation 
results are exactly reproduced. In a forthcoming work we will present results which also 
include exchange diagrams. 

In the present work we do not address the renormalization question by ignoring vacuum 
effects. Here only matter effects are considered up to second order in 6 which includes direct 
as well as exchange graphs. The results are compared with the traditional Hartree and 
Hartree- Fock approximations. Our aim is just to establish the reliability of the method in 
coping with nuclear matter problems. As a byproduct we hope to provide the reader with a 
powerful alternative tool which can be used in investigations aiming to include higher order 
contributions (such as vertex corrections) and vacuum effects. 

Before launching into the actual applications a last remark on how to implement the 
optimized 6 expansion is in order. The standard procedure is to expand the physical quantity 
of interest (P) in orders of 6 starting with the interpolated Lagrangian density C{6). For 
example if P is the energy density (S) one calculates vacuum to vacuum diagrams order by 
order as in perturbation theory using the 0(5°) propagator. In this way 6 labels the diagrams 
contributing to S and improvements will eventually result from the inclusion of higher order 
terms. Alternatively one can obtain an exact expression for S using the energy-momentum 
tensor (T^'^) derived from the original theory C{1). In this case the final expresion for the 
energy density is obtained in terms of the full propagators, which are then evaluated via the 
6 expansion. 

In the present work we adopt the latter prescription for the Walecka model. However, 
the standard prescription is being used in a forthcoming work where the energy density is 
derived perturbatively from the generating functional of the interpolated theory. We will 
then be able to check the equivalence between both prescriptions. 

The paper is organized as follows: in the next section we use the X(j)^ model to pedagog- 
ically introduce the 6 - expansion method. In section 3 we apply the 6- expansion to the 
Walecka model and in section 4 we present our conclusions. 
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2. The model 

To start with, we consider the scalar theory whose Lagrangian density is given by 

£=i(a,0f-imV-^0^ (3) 

To implement the linear S expansion one can consider a general free scalar Lagrangian 
density such as 

>Co = l{d,<Pr - \mlcl? , (4) 

where 

ml = VP? + (5) 

II being an arbitrary mass parameter. Then, according to Eq. (1) one gets 

m = lid,ci>r - \{m' + ii')ct? - 5(^04 - (6) 

The general way the method works becomes clear by looking at the Feynman rules generated 
by C{5). First, the original 0^ vertex has its original Feynman rule —i\ modified to —i5\. 
This minor modification is just a reminder that one is really expanding in orders of the 
artificial parameter 5. Most importantly let us look at the modifications implied by the 
addition of the arbitrary quadratic part. The original bare propagator 

^A(P^) = , , (7) 



becomes at zeroth order in 5 



p2 _ ^2 _ ^2 p2 _ ^2 _|_ 



p2 _ jy^2 _|_ 



(8) 



or 



iA(p^) = — + — -{-i/j,^)— 

pz — 772^ + ^e — rrr + te — + 

+-2 — — — W- + - ' (9) 

pZ _ _|_ pZ _ ^Z _|_ pZ _ ^Z _|_ 

indicating that the term proportional to contained in Cq is entering the theory in a 
nonperturbative way. On the other hand, the piece proportional to S/s^cf)^ is only being 



treated perturbatively as a quadratic vertex (of weight iSfi^). Since only an infinite order 
calculation would be able to compensate for the infinite number of (— i/i^) insertions con- 
tained in Eq. (P) one always ends up with a fj, dependence in any quantity calculated to 
finite order in 6. 

Following the procedure outlined in the Introduction, the final expression for the quantity 
P one wants to evaluate is written in terms of the full propagators which, for the X(f)^ theory, 
is: 

where S(p^) is the self energy. The 6 expansion is then implemented via the substitution: 

tA*ip') ^A\p') = , ^ . , (11) 

where S'^(p^) is calculated perturbatively in powers of S. This implies P = P{iJt) and the 
nonperturbative results are obtained by applying the PMS directly to this quantity, as in 
Eq. (0). 

3. Walecka model 



In this section we consider the Walecka model [0 for nuclear matter. The Lagrangian 
density of the model is given by 

- ^F.^F^'^ + \mlV,V^ , (12) 

where represents the nucleon field operators, and are respectively the field operators 
of the scalar and vector mesons, and F^^ = d^jVu — di^V^. 
We are interested in the energy density of the system: 

g = ^ Jd^x (<^|T°°|^> - <vac|T°°|vac>) , (13) 

where > is the interacting ground-state of nuclear matter, |vac > is the vacuum state 
(zero density), and is the 00 component of energy- momentum tensor T'^'^: 

T^^ = i^Yd^tp + d''(j)d''(f) + d'^VxF^" - g^-'Cw ■ (14) 



Note that we have not used the nucleon equation of motion. 

Next, we express the energy density in terms of full propagators and full self-energies []T3 

^^w = / 1^ Tr {S{k) [7°fc° - irk, - M)] } - ^ / A.(fc) [1 {k^ - ml) - (k' 



(2v 

- <vac|T°°|vac>, 



{e - ml) - [ky 



+ i 



(27 

d^k 



[U,{k)A^{k) + WJ{k)A^,,{k)] 

(15) 



where S{k), Acr{k) and A'^^{k) are respectively the nucleon, scalar- and vector-meson full 
propagators, and the meson self-energies liaik) and 11!^ {k) are given by: 



in,{k) = -ig, / 7^ Tr [S{k + q)V,{p + g, g)^(g)] 



(2vr) 



{2^ f5\k)Uj 7^Tr[5(g)][ , 



(2vr)^ 



(16) 



(fc) = +2^?. / ^ Tr [7'^5(A; + g)C(p + g, g)5(g)] 

- (27r)V(fc) 1^7./ (0iTr[7^5(g)]| |^./ (^Tr[7'^5(g')]} • (17) 

In these equations, the quantities Fj, i = cr, are the full meson- nucleon vertex functions. 
These, in turn, are solutions of Schwinger-Dyson equations that involve higher-order vertex 
functions (or scattering T-matrices). The corresponding bare vertices are given by: 



(18) 
(19) 



It is worth emphasizing that in Eq. (1T5|) direct and exchange contributions as well as 
vertex corrections are included, independent of the order in 5 considered. 

The strategy now is to evaluate the propagators (self-energies) and vertex functions 
according to the perturbative-variational scheme of the optimized 5 expansion discussed in 
the introduction. According to Eq. (|l]), to implement the 5 expansion one needs to introduce 
a £o such that: 



w- 



(20) 



We choose: 
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£o = ^ {il,d^ -Mo)i^ + - ml<p') - \f^.F^^ + \mlV,V^ , (21) 

where 

Mo = M + fi . (22) 

The interpolated Walecka model is then given by: 

£w(5) = Co + 6 [-g^^^^V^ + g„i)(l)^ + /x^'V') • (23) 

Notice that the 5 expansion technique could have also been applied to the meson fields 
explicitly. However, we will eliminate the meson propagators in terms of the nucleon prop- 
agator, using the exact Schwinger- Dyson equations for the meson propagators: 

A,{k) = A^oik) + A,o(A;)n^(A;)A,(A:) , (24) 
Ar(A:) = AZik) + Ai:',{k)U^UkWJ{k) , (25) 



where the meson self-energies are given in Eqs. ([16|)- ([17|) . In this way, meson effects enter 
via the nucleon self-energies. This leaves us with only one unknown parameter, fi, which 
will be fixed by the PMS condition applied to the energy density. As already discussed, the 
implementation of the method will be done via the nucleon propagator, which depends on 
the order in 6 considered. 

Notice also that we could have eliminated the meson-nucleon interaction terms (the terms 
proportional to the meson self-energies in Eq. (^)) by using the exact nucleon Schwinger- 



Dyson equation. This would cancel half of the meson kinetic energies [T^ . In the Appendix, 



we discuss an alternative way to derive the energy density |T^, appropriate to calculations 
up to 0{6^), in which one eliminates from the beginning the meson field operators in favor 
of the nucleon ones. 

At zeroth order in 6, the nucleon self-energy, corresponding to the interpolated La- 
grangian Eq. (|23|) , is obviously zero, i.e., S*^*^) = 0. At this order, the single-particle energy 
is simply given by: 

E{q) = Eo{q) = {f + M^Y' , (26) 
and the nucleon propagator is 



(27) 



with 



q^-M^+ie 



(28) 



S'niq) = il'q, + M,) (f - E{q)) 9 {P, - \q\) , (29) 

the Feynman part of the propagator, Eq. (|28|) corresponding to the vacuum part and Eq. (p9D 
corresponding to the medium part. In what follows we do not consider vacuum contributions. 
At this zeroth-order, meson propagators Acr{k) and ^'^{k) of Eqs. (El)- (12^) are simply: 



A,(A;) = A,o(fc) + A,o(A;)n,(fc)A,o(A;) , 
Ar (A;) = Ko{k) + Kl{m^xAk)Ko{k) 

with the self-energies T\.„{k) and Tl.'^{k) given by: 



n.(fc) = -igi 



2 / d^q 



(27r)4 
(2^ 



Tr 
Tr 



S\k + q)S%q)\+t{27:)H\k){g, 



YS%k + g)7"5°(g) 
d^q 



+ i{27^f5\k) J ^jAYS%q)] g 



d^q' 

(27r)4-L/ - v^yjj i^^-y ^27r)4 
Hence, we obtain for the zeroth order energy density 



TilYS'iq')] 



g{0) ^ ^(0) _^ g{0),dir g{0) ,exc _j_ ^{0),dir _j_ ^{0),exc 



(30) 
(31) 



Tr[5°(g)] , (32) 



(33) 



the following expressions: 



£ 



B — 1 



Pf d^q q^ + MMo 
(27r)3 Eoiq) ' 



1 ^2 

c(0),dir _ __yo_ 

2 ml 



7 



g2 /-Pf 
-7 



rf3^ Mo 

(27r)3Eo(g) 



n 2 



2 ' Jo (27r)3Eo(g) (27r)3Eo(A;) 



-A.([Eo(g)-i5o(A:)]'-(g-fc)' 



X 
X 



- - Ij - [Eo(g) - E,{k)Y/\^{\E,{q) - Eoik)Y - (g - ky 
E^{q)E^{k)-q-k + Ml 



(34) 
(35) 



(36) 
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fiO) 
'-',.1 



2m2 



7 



= gli 



(27r)3 



(37) 



-^^{[E,{q)-E,{k)f-{q-kf) 



(27r)3Eo(g) h (2TifE^{k) 

[E,{q) - Eo(fc)]^A.([Eo(g) - E,{k)f - (g - 
Eo{q)E^{k)-q-k-2Ml\}. 



1 

- - 1 
2 



(38) 



In this expression Aj(A;^), z = a, u; is given by: 



A,(A;2) 



q^ — mj + ie 



(39) 



Note that the term proportional to k^k'^/m'^ in the vector-meson propagator is dropped due 
to the conservation of the baryon current. 

Now we proceed by applying the PMS to S^: 



d£^ 



dS^^ dMo d£^^^ 



0. 



(40) 



d^ dMQ dn dMf) 

At zeroth order in 5, one can see from Eqs. (pO|) and (^) that no interaction between 
mesons and nucleons are considered. Thus, S^^-* = 0. On the other hand, Mq depends 
on /i, vide Eq. (^2]) and it is precisely this parameter, fixed by the PMS condition, which 
introduces all the non perturbative information related to the interactions. Although we are 
working at zeroth order in 5, contributions from direct and exchange terms are included in 
equations (pSD - (^) above. 

Let us first consider the contribution from the direct terms only, which are given by 
Eqs. (^) and (0). Application of the PMS to them yields the following self-consistency 
condition for Mq: 

.2 rPF d^q Mo 



Mn 



mi Jo 



(41) 



(27r)3i?o(g) 

This is exactly the same self-consistency condition for the effective nucleon mass obtained 
by means of the Hartree, or mean-field, approximation. 

Now, application of the PMS to the full energy density leads to a nonlinear equation for 
fi, or equivalently for Mo, which is more complicated than the one of Eq. (^T]). To avoid 
this cumbersome expression, we have chosen to find the minimum of the energy density 
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numerically. In Figure 1 we compare the nucleon binding energy, E /A—M, obtained by using 
only Eq. ( P^ and the direct contributions from Eqs. (^) and (0) (solid line) and coupling 
constants fixed by fitting the binding energy and density of equilibrium nuclear matter, with 
the full binding energy, keeping the same coupling constants (dotted line). The value of the 
coupling constants are g1 = 91.64 and g"^ = 136.2. The masses used in all calculations are 
M = 939 MeV, m„ = 783 MeV and = 550 MeV. We find that the full result coincides 



with those obtained in a relativistic Hartree-Fock calculation p^JT6[1 which we also show 
for comparison (long-dashed line). Note that the dotted and long-dashed lines coincide in 
the figure. Of course, one could renormalize the model parameters to reproduce the bulk 
saturation properties of nuclear matter. This would give us the same couphng constants used 



in the the relativistic Hartree-Fock calculation of Ref. Therefore, the PMS condition on 
the energy density obtained with the zeroth order propagator of the Walecka model is also 
equivalent to the usual Hartree-Fock solution. This is indeed a very interesting result since 
the self-energy expressions are not present and therefore only the exchange contributions to 
the energy density are enough to reproduce, through the minimization of this expression, 
the usual Hartree-Fock result. 

In Figure 2 we compare the results for the effective nucleon mass in nuclear matter as 
a function of Pp obtained from fi. From this figure, it is clear that the results with the 
exchange terms and renormalized constants coincide with the results obtained by using the 
direct terms only. 

Next we check how the previous results change by dressing the nucleon propagator up 
to 0{6^). For this purpose, we start from the calculation of the self-energy. 

For infinite nuclear matter, because of the translational, rotational, parity and time 
reversal invariances, S(g) can be generally written in terms of the unit matrix and the Dirac 
7^ matrices as follows [ p^ : 

S(g)=S^(g)-7^S^(g) 

= S^(gO, |gl) - 7°S°(gO, |g1) + 7 ■ g S^(g°, |g1) ■ (42) 

The self-energy to second-order in delta is given by: 
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• 2 r-2 



S^'\q)K[{p-qf]-i9l5' 



d^q 



j,S^'\q)A^[{p-qfh,, (43) 



(27r)4 "^^^ ' J (27r)4 

where the superscript (2) means second order in 6 and A^- and are given in Eq. (p9|). 
Again, we have made use of the baryon current conservation. 

It is important to notice that this is not a self-consistent equation for although 



its formal similarity with the corresponding Hartree-Fock ones |Tj]. The r.h.s. of this 
equation is expressed in terms of functions calculated at the zeroth-order in 6, as is usual in 
a perturbative calculation. We evaluate Eq. ( P5| ) neglecting the Feynman part of the nucleon 
propagator and considering just ^^(g) given by Eq. ([29|). Because of this, all integrals in 
Eq. (|43|) are finite and can be easily evaluated. The first term in Eq. (|43| ) comes from the 
first order contribution in 6 and must be kept at second order. These expressions are very 
similar to the ones obtained with the Hartree-Fock approximation ||16|. Since there are 
subtle differences, we write them explicitly below. 



Jo [27iy Eq 



+ 



Arr^pJo '"^'^Eoiq) 
gl6^ jPF d^q 





-gyQ.{p,q)-gyQ^{p,q) 



-7- 



1 



(27r) 
dq q 



Att^P Jo 



Pf 



dq q 



-gl5''^a{p,q) + gl5^^^{p,q) 



Air'^p'^ Jo Eo{q) 
where the functions Qi{p, q), $i(p, q),i = cr, ci;, are defined by: 

Ai{p,q) + 2pq 



(44) 
(45) 
(46) 



&i{p,q) = In 



Ai{p,q) - 2pq 



*i(p, q) = '^Mp^ q)Qi{p, g) - 1 



(47) 
(48) 



where 



A{p, q) = f + q' + < - [E{p) - Eo{q)y 



(49) 



One should pay attention to the fact that now the self-energy also carries direct and 
exchange contributions. 

We are in the position to calculate the energy density. We start by defining the following 



auxiliary quantities ||16|| : 
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q =q 



E*iq) = [(r^ + M*\q)]^ , 

and writing the nucleon propagator in the compact form: 
S{q) = Spiq) + Soiq) , 
Spiq) = [1% + M*iq)] ,L/^^^ > 



(50) 



27r 



E*{q) 



SD{q) = [l% + M*{q) 
where E{q) is the single-particle energy: 



5(g°-E(g)h^(P^-|g1) 



(51) 
(52) 

(53) 
(54) 



Note that we have assumed that the nucleon propagator has simple poles with unit residue. 
Within the approximation scheme we are working in this paper, this assumption is satisfied, 
as can be seen below. 

At this order, for the functions Acr{k), A^'^{k), Ilcr{k) and 11^'^ (A;) one has the same 
expressions as in Eqs. (p0|-p3|), where instead of one uses the S above. In what follows, 
vacuum contributions have again been neglected. Hence, we obtain for the energy density 
of nuclear matter the following expression: 



c(2) 



Pf d^q q- q* + MM*{q) 



(27r)3 



E*{q) 



_^ g{2),dir g{2),exc _^ g{2),dir _^ g{2), 



with 



c{2),dir _ _tM£- ^ 
g{2),exc _ ya_^ 



1 gl fPp d^q M*{qj 



(27r)3 E*{q) 
d^q fPF 



(55) 



(56) 



d^k 



X 



2 'Jo {2n)^ E*{q) Jo {2TT)^E*{k) 
1 

- - 1 
2 



1 ^2 

fi2),dir _ 

~ 2 m2 



X 



7 



Pf 



[E{q)-E{k)YA^[{q~ky 
2 



[K[{q-kf] 

q*>'kl + M*{q)M*{k) 



^ d^q 
(27r)= 



(57) 
(58) 



d^k 



(27r)3E*(g) Jo {2'KfE*{k) 



[E{q)-E{k)fA^[{q-kf 



-{A^[{q-kf 



- 2M*{q)M*{k)\] 



(59) 
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These expressions are very similar in form to the ones obtained in the Hartree-Fock approxi- 
mation. However, one should notice that the self-energies are not calculated self-consistently 
as in the Hartree-Fock approximation, rather they are given by Eqs. (|4^) - (^61) , which de- 
pend on Mo, which by its turn, is determined numerically by minimizing the energy density. 
Also, differences are contained in the fermion kinetic energy, the first term in Eq. (|55D, and 
in the factors (^^ — 1^ in Eqs. (|56D and (^). These differences arise because we are not 
using the nucleon Schwinger-Dyson equation. Please refer to the Appendix for an explicit 
derivation in the case one chooses to eliminate the meson field operators from the be- 
ginning. Application of the PMS to the direct contributions present in Eqs. ( p6|) and (^8|), 
calculated only with the direct contributions to the self-energies, yields again the familiar 
Hartree result, i.e., 

M* = M-^,r^^. (60) 

From this result it is straightforward to see that when only the direct terms are considered 
in the energy density and self-energies, the mean-field solution is reproduced at any order in 
5. This result should be compared with the one presented in Ref. where, in the context 
of the effective potential, it was found that the 5- expansion and the expansion are 
identical in the large limit. 

For the full energy density, Eq. (p5|) has to be minimized in terms of /i and this is done 
numerically. This is indeed simpler than the traditional Hartree-Fock procedure, where three 
coupled equations (the self-energy expressions) have to be solved self-consistently. 

We do not present in Figure 1 the 0(5^) binding energy because it would be indistin- 
guishable from the HF one. Instead, for comparison purposes, we present in Table 1 the 
results obtained with the self-consistent Hartree-Fock approximation and the ones with the 
(5° and 5^ expansions. We note that a simple iterative procedure for solving the Hartree-Fock 
equations do not converge for Fermi momenta larger than Pp ~ 1.7 fm~^. Inspection of the 
Table reveals the nice convergence towards the Hartree-Fock approximations of the results 
from 5° to 5^. Moreover, one sees that in order to reproduce the Hartree-Fock results, it is 
enough to use the simple calculation at zeroth order. 

The behavior of M* as a function of the Fermi momentum at this order does not show 
any noticeable difference as compared with the zeroth order results. In Figure 3 we plot 
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the energy density S as a function of fi for Pp = 1.19 fm~^. The sohd hne is obtained 
without the inclusion of the exchange term (the PMS solution in this case is given by 
fi/M = —0.275) and the dashed line gives the full second order density energy (the PMS 
solution is /i/M = —0.35). Recall that if one had an exact solution, the energy density would 
be independent of /i. In this sense it is gratifying to notice that £^ is a very flat function of 
fi. This stability in the value of the energy density as a function of fi is very desirable and 
guarantees that even big changes in the value of /i will not affect physical quantities, as the 
binding energy for instance. 

It is important to point out that although we have obtained the same results for the 
binding energy within the zeroth and second order approximations, this is not true at all 
orders when exchange terms are included. At fourth order in 6, for example, vertex correc- 
tions will appear and the resulting energy density will certainly be different. In this work we 
have opted for neglecting vertex corrections in order to be able to compare our results with 
Hartree and Hartree-Fock results, where they are not included either. If vertex corrections 
are to be included, the full meson-nucleon vertex functions Fj, i = a,u! appearing in Eqs. 



(p!6|) and (|T^ and the full meson propagators will also have to be expanded in orders of 6. 



4. Conclusions 

In the third section of this paper we have utilized the optimized S expansion to study 
medium effects in the Walecka model. We have obtained results quantitatively similar to 
the ones of the usual Hartree-Fock approximation, although the analytical expressions are 
not evidently equivalent. If one neglects the exchange term in the energy density and self- 
energies then clearly the mean-field solution is reproduced at any order. This outcome 
reflects the fact that this perturbative method generates nonperturbative results due, of 
course, of the variational nature of the PMS. 

Figure 1 and Table 1 show that the very simple calculation at zeroth order in 5° already 
provides a very good approximation to the Hartree-Fock results, at least for densities not 
much higher than the normal nuclear matter density. Analytically, this calculation is indeed 
simpler than the usual Hartree-Fock approximation, in view of the perturbative nature of 
the method. Numerically, this calculation is straighforward because no self-consistency has 
to be achieved; one needs only to perform a minimization of the energy with respect to the 
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parameter /x. It is also worth mentioning that, in the Walecka model, the energy density is 
a very flat function of fi and this guarantees that the PMS solution is indeed very stable. 

On the basis of our results, we believe that the optimized S expansion is a very robust 
nonperturbative approximation scheme. Compared with the Hartree-Fock approximation, 
the 6 expansion is very economical because of its perturbative nature. Once the reliability 
of the scheme has been established, one is ready to proceed to other interesting applications. 
These include vertex and, obviously, vacuum effects that include exchange corrections |]T3[ . 
In view of our results, we can proceed by including vertex corrections in the energy density 
and still maintaining the nucleon propagator at zeroth order in 6. The study of the vacuum in 
the Walecka model is an important issue since one needs to know the limits of applicability 
of such model to high densities and/or temperatures before quark and gluon degrees of 
freedom have to be invoked. Particularly interesting is the renormalization of exchange 
diagrams which should become simplified in the present approach as compared with the 
Hartree-Fock scheme |]T^, since at each order in 6, only a finite number of diagrams has to 
be taken into account. 
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TABLES 

TABLE I. Results for E/A - M (MeV) as a function of Pp (fm"^) calculated with the 
Hartree-Fock (Ehf) approximation and with the 6^ and S'^ {Ep) expansions. 



p 

-TF 


TP 


TP „ 




U.Uo 


u.uzyiuyo 


u.uzyiuyo 


u.uzyiuyo 


f) 25 


5S1 75S0 


5S1 752Q 


531 7530 


0.50 


1.1637572 


1.1637562 


1.1637569 


0.75 


0.5920830 


0.5920533 


0.5920831 


1.00 


-1.8171553 


-1.8175904 


-1.8171548 


1.25 


-4.0899150 


-4.0937442 


-4.0898972 


1.50 


5.7406683 


5.7109528 


5.7409897 


1.65 


31.0297464 


30.9322795 


31.0312656 
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FIGURES 

FIG. 1. Pp dependence of the binding energy of the Walecka model at zeroth order in 5. The 
solid line represents the direct contributions only (Eqs.(34),(35),(37)). The dotted and long-dashed 
lines give the full binding energy and the Hartree-Fock solution respectively, both determined with 
the same coupling constants used in the solid line solution (note that the lines are coincident). 

FIG. 2. Zeroth order nucleon effective mass Mq as a function of Pp- The solid curve is the 
result obtained without the exchange term and the dashed curve is the result using the full energy 
density. 

FIG. 3. jj, dependence of the energy density for the Walecka model at second order in 5, calcu- 
lated at Pp = 1.19 fm~^. The solid line gives the solution when the exchange term is not included. 
The dashed line gives the full solution. 
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Appendix 



For completeness, in this Appendix we employ the method used in Ref. [0 for obtaining 
the general expression for the energy density in terms of the nucleon propagator, valid up 
to order 0{6'^). 



The energy-momentum tensor is defined by Eq. ( p!2D with Cw defined in Eq. (|Tj): 

T{^- = + Tr + (61) 

with 

T^" = i^Yd^i^ - g^^^ihad" - M)^ , (62) 
Tr = - g'-'' [^g^n^ - \ml<p' + \d^(t^d^<p) , (63) 

Tr = d'VxF^^^ - g'^'' (-(7.7a V^" + - ^F^pF'^^^ . (64) 

Let us concentrate on T^'^ . The Euler-Lagrange equation for the scalar-meson field 
equation is: 

(a^a^ + ml)<f) = gM ■ (65) 
This equation can formally be integrated as: 

(j)[x) = 0°(x) - ga J d'^yKix - y)i'{y)^{y) , (66) 

where 0° is the solution of the homegeneous equation and A^{x) is given by: 

A.(x) = / TSl-^^-^e-- = / ^ A.(g^) e"-. (67) 



(27r)4 g2 _ ^2 + J (27r)4 

From Eq. (^) (note that the first term (j)^{x) does not contribute) we have: 



ga < >= -g^ J d yK{x - y) < i)a{x)i)i3{y)'4)p{y)i)o,{x) > . (68) 
With the help of Wick's contraction technique we get 



/ (0^""^'^'"'^ / (0i^"^'^'"^^Tr[5(g)5(A:)] 
21 



(69) 



that can finally be written as 



9l 



g„ < ^(f)ilj >= ^ 



(2vr) 



■Tr[5(g)] 



9l 



{2txY (27r)^ 



Ti[S{q + k)S{k)]K{q^) . (70) 



The third term in Eq. ( pB| ) above can be written as 

1 1 r 

- -ml < 4>^ >= --mlgl j d'^yd^zA^ix - y)A,{x - z) < ^Pa{y)ijo,iy)i'piz)^f3{z) > . (71) 

Using again Wick's technique 
1 



- -ml < (h^ > 



\mlgl j rf^yrfS^^e--(-^)A.(/)e--(^-^)A.(g^ 
d^k' 



d'^k 
d^k 



IY[5(fc)] 

-ik(z-y) 



(2vr)4 
d^'k' 



Tr[S{k')] 



(2vr) 



e-i'''(y-')TT[S{k)S{k')] 



(72) 



and finally 



-ml < d)'^ > 



2 ml 
m„q„ 

2 (T J tJ 



dS d^k 



AM^)i:i[S{q + k)S{k)]K{q^) 



(73) 



(27r)4 (27r)4 

Following the same procedure for d^(f)d^(j) and du<f)d'^(f), one obtains after adding all terms: 



-lOO 



+ 



l9l 
2 ml 



dS 
(2^ 



Tr^(g) 



■ ' 2 [ d\ 
+ (2^ 



d'^q d'^k 



(27r)4 (27r) 



Ti[S{q + k)S{k)]A„{q^ 



X 



(g°)^A.(g^) + ^(g^ + m^)A.(g^) 



(74) 



The same can be repeated for the vector-meson field. This shows the equivalence up 
to 0{5'^) of the perturbative solution of the Schwinger-Dyson equations and the method of 
elimination of the meson field equations from the beginning 



22 



1.0 



0.5 



0.0 



